Magnetic and thermodynamic properties of cobalt doped iron pyrite: Griffiths Phase 

in a magnetic semiconductor. 
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Doping of the band insulator FeS2 with Co on the Fe site introduces a small density of itin- 
erant carriers and magnetic moments. The lattice constant, AC and DC magnetic susceptibility, 
magnetization, and specific heat have been measured over the < a:: < 0.085 range of Co concen- 
tration. The variation of the AC susceptibility with hydrostatic pressure has also been measured in 
a small number of our samples. All of these quantities show systematic variation with x including 
a paramagnetic to disordered ferromagnetic transition ai x — 0.007 ± 0.002. A detailed analysis 
of the changes with temperature and magnetic field reveal small power law dependencies at low 
temperatures for samples near the critical concentration for magnetism, and just above the Curie 
temperature at higher x. In addition, the magnetic susceptibility and specific heat are non-analytic 
around H = displaying an extraordinarily sharp field dependence in this same temperature range. 
We interpret this behavior as due to the formation of Griffiths phases that result from the quenched 
disorder inherent in a doped semiconductor. 
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PACS numbers: 75.50.Pp, 75.40.-s, 75.20.Hr 



I. INTRODUCTION AND MOTIVATION 

The desire for ever increasing computational speeds 
has led to the conception of spintronics technologies that 
make use of both the charge and spin properties of elec- 
trical charge carriers in solidsii^. Reahzation of these 
devices for useful technologies requires the discovery and 
development of materials that have easily controllable 
electronic and magnetic properties to the same extent 
that semiconductors allow control over charge proper- 
ties. This demand has led to an enormous interest 
in magnetic semiconductors to investigate if traditional 
semiconductors can be manipulated by chemical substi- 
tution of transition metal elements, or other creative 
means, into offering control over the spin properties of 
electrons^'^'^'i'^. Because so many of the suggested de- 
vices rely on the production and detection of spin polar- 
ized currents, magnetic semiconductors vi^ith Curie tem- 
peratures above room temperature that are compatible 
with silicon technology are highly desired. However, this 
goal has thus far not been met and it has become in- 
creasingly apparent that to make progress, fundamental 
knowledge and understanding of the electronic and mag- 
netic properties of magnetic semiconductors are neces- 
sary. Although some headway in this direction has been 
made, the complexity of disordered materials with strong 
interactions makes magnetic semiconductors difficult to 
model and understand^ d^d^ . 

In a recent paper we presented the results of an explo- 
ration into the magnetic and electronic properties of one 
such magnetic semiconducting compound, Fei_a;Coa;S2 
based on the diamagnetic insulating parent compound 
iron pyrite (FeS2)i^. We found that the insulator-to- 
metal transition that occurs for low concentrations of Co 
doping was followed by a transition from a paramagnet to 



a disordered ferromagnet first apparent for small dopant 
concentrations, x, at very low temperatures. Increased 
Co doping led to an increased transition temperature and 
to a more ordered ferromagnetic phasei^^'^^'^ii^. The 
main conclusion of our paper was that the quenched dis- 
order associated with the chemical substitution has im- 
portant consequences on the formation of the magnetic 
state at small x and on the finite temperature phase tran- 
sitions at larger x. Unusual power law temperature, T, 
and magnetic field, H, dependencies of the magnetic sus- 
ceptibility, magnetization, and specific heat were mea- 
sured at low-T for x near the critical concentration for 
nucleation of a magnetic phase, well as at T's just 

above the Curie temperature, Tc for x > Xc- These were 
attributed to large, rare disorder fluctuations consisting 
of regions of local magnetic order which fluctuate as a 
single large magnetic moment and that grow in size as 
T ^ or T ^ Tc. Fluctuating rare region phenom- 
ena associated with phase transitions in systems with 
quenched disorder are known under the heading of Grif- 
fiths phasesi^J^i^o. 

Theorists have found the subject of Griffiths phases 
fascinating since Griffiths discovered the non-intuitive 
feature in his calculations that, under certain condi- 
tions, exponentially rare regions in disordered systems 
can dominate the thermodynamic response when in prox- 
imity to a phase transitioni*'^*'. This is in contrast to 
many aspects of phase transitions where disorder can be 
treated perturbatively or simply ignored. These ideas 
have drawn even more interest as evidence for non- Fermi 
liquid behavior in disordered systems has been discov- 
ered while explanations remain tentative^!. While it is 
clear that the disorder and strong electron correlations 
are responsible for the breakdown of Fermi liquid behav- 
ior in some cases, the microscopic origins remain poorly 
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understoodiS,. In this context, the unusual power law be- 
havior associated with the formation of Griffiths phases 
have been discussed as one mechanism for disorder driven 
non-Fermi liquid behavio r^^i^^ . Experimental evidence 
that Griffiths phases are observable in condensed mat- 
ter systems has only recently been reported for a small 
number of systems2ii2^i2£i21i2^, and in our previous publi- 
cation we presented evidence for the importance of Grif- 
fiths phase formation near the critical concentration for 
magnetic ordering in the magnetic semiconducting sys- 
tem Fei_a;C0a;S2. The purpose of the present paper is 
to present a more thorough rendering of our magnetic 
susceptibility, magnetization, and specific heat experi- 
ments and analysis, as well as to present the results of 
some more recent AC magnetic susceptibility measure- 
ments made at ambient and applied pressures of up to 7 
kbar. In a separate accompanying article we present the 
results of our Hall effect and resistivity measurements of 
these same materials^^. 

We chose to investigate Co doped FeS2 since FeS2 is a 
simple diamagnetic band insulator with the same pyrite 
crystal structure as C0S2 , a itinerant ferromagneti^illi^. 
The two compounds had been shown earlier to be misci- 
ble such as to form the substitutional series Fei_a;Coa;S2 
without observation of second phases or any phase seg- 
regation that occurs in many magnetic semiconducting 
compounds. As such, it allowed the investigation of the 
transition from diamagnetic insulator to a paramagnetic 
metal, then to a ferromagnetic metal with Co substitu- 
tion. In addition, band structure calculations have pre- 
dicted that at large x, 0.5 < a; < 0.95, Fei_a;C0;ES2 
is half metallic having a fully spin polarized ground 
statei^i^^i^ 

We begin our discussion by presenting the details of 
the experimental techniques that we have used to inves- 
tigate the extrinsic properties of Fei_a;Coa;S2. This is 
followed by a presentation of the magnetic susceptibil- 
ity, magnetization, and specific heat over a wide range 
of X covering insulating, paramagnetic metallic, and fer- 
romagnetic metallic behaviors. In section HVl we explore 
in detail the behavior of our crystals nearest at low 
temperatures where the power-law in T and H behaviors 
are found. This is followed by a discussion of the mag- 
netic properties of crystals with x > Xc in proximity to 
Tc displaying similar behavior. We conclude the paper 
with a summary and discussion of our results. 



II. EXPERIMENTAL DETAILS 

Single crystals of Fei-xCo^Si were synthesized from 
high purity starting materials, including Fe powder (Alfa 
Aesar 99.998%), Co powder (Alfa Aesar 99.998%), and 
sulfur (Alfa Aesar 99.999%). These materials were sealed 
in a 16 mm diameter quartz tube under vacuum. Iodine 
(Alfa Aesar 99.99-1-%) was added to the tube as the trans- 
port agent at a concentration level of 2 to 5 gm/cm^. 
The tubes were then placed in a small, single zone hor- 



izontal tube furnace (Lindberg/Blue Model 55035) and 
heated to a constant temperature. Insulating plugs were 
removed from the ends of the furnace to create a natural 
temperature gradient of ^ 50 °C. The charge was placed 
at the center of the furnace (hot zone), and single crys- 
tals would begin to form at the cold end of the tube in 
approximately 5 days, with a maximum yield by no later 
than 14 days. Good crystals were produced for tempera- 
tures in the range of 650 to 850 °C. Crystals were etched 
in HCl to remove any remaining flux. 

Initial characterization included single crystal X-ray 
diffraction. Single crystal fragments were glued to a 
glass fiber and mounted on the goniometer of a Nonius 
Kappa CCD diffractometer equipped with Mo-Kq, radi- 
ation (wavelength — 0.71073A). All X-ray structure de- 
terminations are consistent with a Pyrite crystal struc- 
ture with no anomalous intensities that could be inter- 
preted as evidence for a second phase. The model of the 
structure was refined to a level consistent with a fully 
solved crystal structure. We observe no evidence for dis- 
order, including that associated with sulfur deficiencies, 
in the crystal structure of all our crystals. The lattice 
constants, a, shown in Fig. [1] were determined by the 
X-ray diffraction measurements. As displayed in the fig- 
ure, a increases systematically with Co concentration, 
X in Fei_2:Co2;S2, beyond the lattice constant for pure 
FeS2, a = 0.54165 nm, and is comparable to the mea- 
surements of Refs.[llllll|I3. The increase in lattice 
constant with x is consistent with Vegard's law and the 
idea that Co replaces Fe within the pyrite crystal struc- 
ture. Energy-dispersive X-ray microanalysis (EDX) on 
a JEOL scanning electron microscope equipped with a 
Kevex Si(Li) detector was performed to check the stoi- 
chiometry of our samples. These revealed that the Co 
concentration was about 0.7 ± 0.1 times the nominal Co 
concentration of the starting materials, while the sum of 
the Co and Fe densities was within error, ±2%, of one 
half the sulfur density. 

Magnetization, M, and magnetic susceptibility (both 
AC and DC) measurements were performed in a Quan- 
tum Design superconducting quantum interference de- 
vice (SQUID) magnetometer (MPMS) for temperatures, 
T > 1.8 K and DC magnetic fields, H, of up to 5 T. 
AC susceptibility measurements were carried out over a 
range of excitation fields from 0.1 to 1.5 Oe and excita- 
tion frequencies from 1 Hz to 1 kHz. In addition, AC 
magnetic susceptibility measurements were performed in 
a dilution refrigerator above 50 mK at = 0, a frequency 
of 1 kHz, and an excitation field of 1.5 Oe. The satura- 
tion magnetization was determined by DC magnetization 
measurements at 1.8 K and high fields (up to 5 T). These 
measurements were consistent with our EDX measure- 
ments showing a magnetic moment density of 0.7 ± 0.1 
times the nominal Co concentration. Because M mea- 
surements are easily performed and highly reproducible, 
we have subsequently used the saturated magnetization 
to determine the Co concentration of our samples. Thus, 
throughout this manuscript the stoichiometry of the sam- 
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FIG. 1: (Color online) Lattice constants. Lattice con- 
stant, a, for a wide range of cobalt concentration, x, de- 
termined by single crystal X-ray diffraction for our samples 
(labeled Guo) along with comparisons with those previously 
determine d^^'^^'^^'^^ . Black line is a linear x dependence of a 
from FeS2 through to C0S2 . 



pies noted in the figures and text was determined in this 
manner. The variations in the saturated magnetization 
for crystals from the same growth batch was measured 
to be ±10% of the average value. Therefore, we report 
X determined from measurements of M at high field for 
the crystals used in each of our measurements. Wherever 
possible, the same crystal was employed for several differ- 
ent measurement types, particularly for the measurement 
of the AC susceptibility and specific heat below 1 K. 

The magnetic susceptibility of several of our crystals 
was measured with applied pressure ranging from am- 
bient to 7 kbar in a beryllium-copper pressure cell^. 
Our pressure cell was based on previously published 
designs^i24i35_ r^-^e pressure cell was machined from 
beryllium-copper alloy 25 since this material has a low 
magnetic susceptibility and superior mechanical proper- 
ties. After machining, the cell was annealed in air at 
a temperature of 335 °C for 2 hours. We designed our 
pressure cell to have a length of 12 cm to reduce the back- 
ground seen by the SQUID pickup coils originating from 
the ends of the cell. The pistons are made from high pu- 
rity Technox 3000 zirconia rods from Dynamic Ceramic. 
This resulted in a small, diamagnetic, (~ —7 x 10~^ emu 
at 50 Oe field) background signal due to the lack of zir- 
conia at the sample position. This background depended 
slightly on the distance between the ends of the two zir- 
conia rods which, in turn, depended on the sample thick- 
ness, increasing for larger separations. The sample was 
placed inside a Teflon container which was sealed with a 
beryllium-copper cap and retainer. The Teflon container 
is filled with silicon oil (Dow Corning 704) as the pressure 
transmission medium (see for example Ref. |3^ for more 
on silicon oil as a pressure transmission medium). The 
pressure at the sample position was determined by plac- 
ing a small (0.5 mg) (Vq. 99 Tio. 01)203 crystal in the Teflon 
container for use as a manometer'^^. This material makes 



an excellent manometer as it displays a readily measur- 
able metal-to-insulator transition which varies linearly in 
T with pressure for pressures of up to 10 kbar—. 

The specific heat was measured using a standard ther- 
mal relaxation method in a Quantum Design Physical 
Property Measurement System (PPMS) above 2 K and 
a dilution refrigerator above 100 mK. For measurements 
below 2 K single crystals were polished fiat and mounted 
on a Sapphire substrate with GE varnish. The oppo- 
site side of the 0.25 mm thick Sapphire substrate con- 
tained a RuO resistance thermometer which was notched 
to provide a thermometer-heater pair in good thermal 
contact with the sample. The thermometer was cycled 
from 300 K to 4.2 K many times to ensure repeatability of 
the resistance at low temperatures'^^. The thermometer 
was calibrated by comparison to a Lakeshore calibrated 
germanium resistance thermometer. Four resistive wires 
provided electrical and thermal contact of the thermome- 
ter and heater to the mixing chamber of the dilution re- 
frigerator. We carefully checked at each measurement 
temperature that the time constant for thermalization of 
all elements on the Sapphire substrate was much smaller 
than the thermal time constant for thermalization to the 
cryostat. Specific heat measurements were carried out in 
magnetic fields between and 6 T with all thermome- 
ters calibrated at each magnetic field. The specific heat 
of the addenda was measured separately and carefully 
subtracted from the data. 



III. OVERVIEW OF EXPERIMENTAL RESULTS 

A. DC susceptibility and Magnetization 
Measurements 

In order to establish the magnetic ground states and 
critical temperatures of the Fei_2;Coa;S2 chemical sub- 
stitution series we have measured the DC magnetic sus- 
ceptibility, xdCi and magnetization, M, of our crystals. 
In Fig. [2] we plot xdc of a dozen representative samples 
with X ranging from nominally pure FeS2 to x — 0.055. 
Our X = crystal has a diamagnetic susceptibility that 
changes little as a function of temperature. The small 
size of the impurity related Curie- Weiss contribution to 
the susceptibility of this crystal observable below 20 K, 
an increase of 10^^ from 20 K down to 2 K, is consistent 
with an impurity concentration of 3 x 10^^ cm~^. Not 
only does this establish the quality of our pyrite crystals, 
it is compelling because it reaffirms that the parent com- 
pound is non-magnetic despite being 1/3 iron^^'^**. The 
low temperature xdc systematically increases with x by 
over 4 orders of magnitude for our a; = 3 x 10~* to our 
X — 0.055 samples. At the same time, the T-dependence 
changes from being almost Curie- like at low- a;, to display- 
ing a broad, almost constant region of xdc below 5 K for 
X > 0.035. At higher temperatures xdc is less system- 
atic in X with some crystals displaying almost constant 
Xdc above 100 K suggestive of a significant Pauli sus- 
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ceptibility, while others show a more T-dependent xdc 
to higher temperatures crossing through xdc of samples 
with smaller x. Because the Pauli susceptibility is pro- 
portional to the conduction electron density of states, we 
interpret the variation of our measured xdc above 100 K 
as revealing a carrier density and, perhaps, a mass that 
is highly sensitive to the disorder. Because of this change 
from a systematic low-temperature xdc to a less system- 
atic high-T Xdc, we conclude that in the two tempera- 
ture ranges Xdc is dominated by two different electron 
populations; a more local population of electron moments 
dominates at low-T while the more itinerant electrons 
dominate xdc at T > 100 K. 



in the data, we have quantified the Curie Weiss be- 
havior by fitting a linear T-dependence to the data as 
shown in the figure. We interpret the intercept of the 
fits with the T-axis as Qyp and the slope as 1/CC 
where CC — n{gijLB)^J{J+l)/'ikB with n the density of 
magnetic moments of size J, g the gyromagnetic factor, 
ks Boltzmann's constant, and the Bohr magneton. 
Both Q^F and CC are observed to increase systemati- 
cally with a; as can be seen in Fig. [4] where we plot the 
best fit values of the quantity ng'^J{J+l) = 3CCkB/ 
and ■ It is also interesting to note that for small 
X, X < 0.01, '^^P is less than zero revealing a small 
antiferromagnetic interaction of the moments for small 
Co doping. This is perhaps not expected since at only 
slightly larger Co concentrations a disordered ferromag- 
netic ground state results. 



FIG. 2: (Color online) DC magnetic susceptibility. Tem- 
perature, T, and cobalt concentration, x, dependence of the 
DC magnetic susceptibility for several representative crystals. 
Note that the data for the a; = 0, nominally pure sample, 
is diamagnetic and thus has been multiplied by —1 for dis- 
play purposes. Magnetic fields applied include 10 kOe for the 
a; = 3 X 10~*, 7 x 10"", and 0.004 crystals; 2 kOe for the 
a; = 0, 0.012, and 0.035 crystals; 1 kOe for the x = 0.005, 
0.007, and 0.045 crystals; and 50 Oe for the x = 0.031, 0.05, 
and 0.055 crystals. Measurement fields were chosen so that 
noise levels were acceptable given the size of the susceptibility 
signal for each individual crystal. 

To gain insight into the character of the magnetic mo- 
ments responsible for the low-T upturn in xdc we have 
plotted in Fig. [3] the inverse DC magnetic susceptibility 
after subtraction of the temperature independent Pauli 
susceptibility, xp, ^/{xdc ~ Xp) a function of T. 
This is standard procedure for comparing the suscepti- 
bihty data to a Curie- Weiss form, x = CC/{T - O^^) 
where CC is the Curie constant and Q^P is the Weiss 
temperature. Although some nonlinearity is apparent 



FIG. 3: (Color online) Inverse DC magnetic susceptibility, a) 
and b) Temperature, T, and cobalt concentration, x, depen- 
dence of the inverse DC magnetic susceptibility, 1/xdc, after 
subtraction of the temperature independent contribution, for 
several representative samples at the same magnetic fields as 
in Fig. [5] The constant susceptibility is attributed to the 
Pauli susceptibility, Xp our crystals. Inset to frame a) dis- 
plays 1/{xdc — Xp) below 100 K for our x = 3 x 10~^ and 
X = 7 X 10~^ where the signal to noise level is such that Curie 
behavior is apparent in an inverse susceptibility plot. Error 
bars are plotted only for every 5th (3rd in inset) data point 
for display purposes. 

We note that the Curie constant increases with x be- 
yond what would be expected for paramagnetic spin-1/2 
magnetic moments for all samples measured. If we as- 
sume that the density of magnetic moments is equal to 
the Co density, and that g = 2 we would conclude that 
the magnetic moment increases from just above J = 0.5 
to between J = 1 and J = 2 for x > 0.03. This is the 
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FIG. 4: (Color online) Curie constant, Curie temperature 
and Weiss temperature, a) The Curie constant, ng^J{J + 
1), determined from fits of a Curie- Weiss form added to 
a temperature independent Pauli contribution, xpi X ~ 
n{giiB)^ JJ{+l)/3kBT + xp, to the DC and AC magnetic 
susceptibility. Lines indicate values of the Curie constant for 
various J's when n is set equal to the cobalt concentration of 
the crystals and g is set equal to 2. b) The Curie temperature, 
Tc as defined by the temperature of the maximum in the real 
part of the AC susceptibility, max x'j and through a mean 
field Arrott analysis, Tc Arrott (see text for details). The 
temperature of the maximum in the imaginary part of the 
AC susceptibility, max x" is included in the figure for com- 
parison to max x' ■ Weiss temperatures as determined from 
the temperature dependence of the AC, O^, and DC, Qw^ , 
susceptibility are also plotted. Lines are fits of a linear max 
x'{x) and <d^' (x) dependence to the data. 



first indication that the magnetic moments likely form 
strongly interacting clusters over the temperature scale 
where the Curie constants were determined. The mea- 
sured Curie constant, in this case, most likely reflects the 
fluctuations of a smaller density of much larger composite 
magnetic moments. 

The magnetic properties of our samples was further 
explored by measuring M at 1.8 K in flelds of up to 5 T. 
Fig. [5^ displays the enormous increase in M with x for 
a range of X between 3 x lO""* and 0.085. In Fig. [5]b 
we display M, in fiB per Co dopant, to demonstrate the 
normalization used to determine the Co density of our 
crystals. The form of the magnetization can be seen to 
evolve from that of a Brillouin function, Bj, at small x 
to that of a ferromagnetically ordered material at larger 

X. 



FIG. 5: (Color online) Magnetization, a) The magnetization 
vs external magnetic field // at 1.8 K for several of our crystals 
with Co concentrations identified in frame b. b) The magne- 
tization normalized so that the high field values approach 1. 
This normalization was used along with the EDX data to de- 
termine X. The red line is a Brillouin function for J = 1/2 
and a gyromagnetic ratio of 2. 



The temperature and fleld dependence of M, such as 
that shown in Fig. [51 is commonly compared to a mean 
fleld form to determine Tc of magnetic materials. In 
Fig. [6] we display this comparison of AI {H, T) for two 
representative samples. The mean field form for M is 
typically written as M{H) = ngjfis Bj{gJfj,B [H + 
r]M)/kBT) with 77 the constant parameterizing the 
strength of the molecular mean field. From this start- 
ing point an expansion for H/M can be written as 
H/M = (1/xo) + aM"^ + bM'^ + ... where xo is the initial 
susceptibility (x(i? = 0)). Thus, a plot of vs. H/M 
should be a straight line for a range of fields such that the 
higher terms can be ignored'^ This is commonly referred 
to as an Arrott plot and is demonstrated by Fig.[n]frames 
a and b showing the results of such an analysis for two 
samples near the critical Co concentration for ferromag- 
netism. It is clear that the region of linearity of in 
H/M is limited, particularly at low temperatures. Frame 
c of the figure displays a compilation of 1 /xo values deter- 
mined from the linear fits, such as those in frames a and 
b. We have plotted the values of 1/xo vs. T^ since the 
Stoner-Wohlfarth theory of itinerant magnetism predicts 
a T^ dependence of 1/xo near the transition tempera- 
ture. Linear fits to the low-T data are used to determine 
the T^ intercept and this is interpreted as the mean-field 
Tc- The values for Tc determined in this way are plotted 
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in Fig.|4]b for comparison to ■ Modified Arrott anal- 
ysis's having different dependencies of H/M on M were 
also performed without significant changes to the criti- 
cal temperatures determined from the above analysis^. 
The conclusion from this mean-field analysis is that the 
Curie temperature is less than 1/20^*-^ for these samples 
reflecting the importance of disorder on the formation of 
the magnetic state and that there is a critical Co con- 
centration for ferromagnetism near x = 0.03. However, a 
simple mean field analysis is likely to be inaccurate when 
Tc is close to zero. 



FIG. 6: (Color online) Arrott plots. Standard mean field 
analysis for determining the Curie temperature, Tc, of ferro- 
magnetic materials. In a) and b) we plot the square of the 
magnetization, M, vs. the external field, H , divided by the M. 
Tc is the temperature where the extrapolation of the data to 
H = goes through 0. a) Arrott plot for an a; = 0.031 crys- 
tal. Data are at temperatures of 1.8 (red filled diamonds), 
2.0 (orange squares), 2.25 (yellow bullets), 2.5 (yellow-green 
open triangles), 2.75 (green filled squares), 3.0 (green open di- 
amonds), 3.25 (light blue filled triangles), 3.5 (blue asterisks), 
3.75 (blue- violet circles), 4.0 (purple pluses), 4.25 (violet open 
right pointing triangles), 4.5 (dark violet X's), 4.75 (black 
left open left pointing triangles), 5.0 (red open squares), 5.5 
(orange filled diamonds), 6.0 (yellow open triangles), 6.5 (yel- 
low green bullets), 7.0 (green open diamonds), 7.5 (blue-green 
filled squares), 8.0 (light blue asterisks), 8.5 (blue filled tri- 
angles), 9.0 (blue pluses), 9.5 (blue- violet circles), and 10 K 
(violet X's). Dashed lines are linear fits to data between 1 
and 15 kOe. b) Arrott plot for an x =0.031 crystal. Tem- 
peratures and symbols same as in frame a. Dashed lines are 
linear fits to data between 1 and 20 kOe. c) Plot of H/M axis 
intercepts, Xo^y '^^ linear fits to data in frames a and b vs 
for several crystals identified in the figure. The intercepts 
of these data with the axis is interpreted as the mean-field 
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B. AC Susceptibility Measurements 

Although the DC susceptibility and magnetization give 
us a good initial indication of the formation of a ferro- 
magnetic phase with the substitution of Co for Fe in FeS2, 
the very subtle changes that occur near Xc are best ex- 
plored at very low magnetic fields where good quality DC 
magnetization signals are difficult to obtain. Therefore, 
we have made very careful AC susceptibility measure- 
ments of our crystals using excitation fields of 1.5 Oe or 
less as shown in Fig.[71 In agreement with our DC results, 
the AC susceptibility shows a dramatic increase with x. 
However, there are two important differences between our 
DC and AC susceptibility apparent in the figures. First, 
the real part of the AC susceptibility, x', is noticeably 
larger than xdc for all x > crystals measured. Sec- 
ond, distinct maxima appear in both x' and, at slightly 
higher T's, the imaginary part part of the AC suscep- 
tibility, x", as can be seen in frames a and b. As is 
common practice, we preliminarily interpret the temper- 
ature of the maxima in x' as the critical temperature, 
Tc, for the development of a magnetic, ferromagnetic or 
spin glass, phase. In Fig. we plot Tc determined in 
this manner (max x') along with the temperature of the 
maximum in x" (max x") for comparison to our analy- 
sis of Xdc and M{H). In contrast to our xdc results, 
our very low-T measurements (down to 50 mK) reveal a 
peak in x' for all crystals with x > 0.007 indicating that 
Xc — 0.007 ± 0.002. The contrasting results in our AC 
and DC magnetic susceptibility indicate that the mag- 
netic state of these crystals are extraordinarily sensitive 
to magnetic field as we explicitly demonstrate in section 
El below. 

In frame c of Fig. [7] is plotted against temperature 
to display the Curie- Weiss-like behavior of x' ■ Here, the 
general Curie- Weiss-like trend is displayed although sig- 
nificant nonlinearity and structure are also apparent in 
these data. Linear fits of these data were used to deter- 
mine the Weiss temperatures, &yi/^ , and Curie constants 
displayed in Fig. 01 There appears to be good compati- 
bility of both Qyf^ with the values determined from the 
DC values, however, the Curie constant appears gener- 
ally smaller for our AC analysis than for our DC val- 
ues for X > 0.02. The Curie constants determined from 
both methods indicate a fluctuating magnetic moment of 
J > 1/2 per Co dopant indicative of magnetic cluster for- 
mation. Despite the general linear behavior of 1/x' well 
above Tc, there is a noticeable departure from this behav- 
ior at T < 48^ that can best be seen for our x = 0.024 
sample in frame c of Fig. [71 Fits of l/x'(T) for T > 
15 K indicate = 3 K, however, at T < 15 K 1/x' 

tends toward much smaller values. This tendency for x' 
to increase at a rate faster than the Curie- Weiss behavior 
indicates local ferromagnetism and further supports our 
observation that magnetic clusters form at temperatures 
above T^. A more complete description of this behavior 
is presented in section |Vl of this paper. 

The effect of hydrostatic pressure, P, on the magneti- 



FIG. 7: (Color online) AC susceptibility, a) The real part 
of the AC magnetic susceptibility, x' of several of our crys- 
tals with Co concentrations identified in frame b. The AC 
excitation field was 1.5 Oe and the measurement frequency 
was 1.0 kHz for the x =0.005, 0.007, 0.012, 0.024, and 0.045 
crystals, 1.0 Oe and 10 Hz for the a;=0.029 crystal, and 1.0 
Oe and 100 Hz for the 1=0.05 crystal, b) Imaginary part of 
the AC susceptibility, x"i for the same crystals in frame a. 
Co concentrations are identified in the figure, c) 1/x' demon- 
strating Curie Weiss-like behavior. Data and symbols same 
as in frame a and b. 



zation of our crystals is demonstrated in Fig. [51 While 
the magnitude of x' is depressed by the application of 
pressures of order of a kbar, there is little change to the 
temperature of the peak in x' ■ The changes that occur 
with P are demonstrated in frame b of the figure where 
the temperature of the peak in x', taken as our definition 
of Tc, is plotted along with the results of fitting a Curie 
Weiss form to x'{T) above 15 K. Both the Curie constant 
and the Weiss temperature are reduced as P is applied 
showing that P reduces both the average size of mag- 
netic clusters as well as the average interaction energies 
between magnetic moments evident in the P = mea- 
surements. The Curie constant shown in Fig. |8lis consis- 
tent with a reduction of the density of J = 1/2 magnetic 
moments from 3 x 10^^ to 1.7 x 10^^ cm^'^ as the pressure 
is increased from ambient to 6 kbar. Equivalently we can 
assume that the density of moments remains equal to the 
Co density and that the changes in the Curie constant 
reflect a reduction of J from 1.3 to 0.9 with P. Our 
data could also be interpreted as a general trend toward 
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a stronger Kondo coupling with P as has been shown to 
be typical in the case of f-elcctron compounds^. 



contribution to the specific heat at temperatures well 
below the Debye temperature. While /? = 2.17 x 10~^ 
J/mol describes the phonon contribution of all our 
samples well, there are significant changes to the elec- 
tronic contribution to C'(T) with x. It is clear that 7 of 
the nominally pure FeS2 crystal is consistent with zero 
and that 7 increases with x. However, there are contri- 
butions to the specific heat that are not described by the 
standard form. For example, at the lowest temperatures 
an upturn in C(T)/T with decreasing T is apparent in 
all of the Co doped samples. Several of the data sets also 
show an increased contribution between 10 and 70 K that 
is not described by a simple constant. It is most likely 
that the differences from the standard metallic form are 
due to the formation of the magnetic ground state with 
Co doping that was apparent in the magnetic suscepti- 
bility. 



FIG. 8: (Color online) Pressure dependence of the AC mag- 
netic susceptibility, a) The pressure, P and temperature, T, 
dependence of the real part of the AC magnetic susceptibil- 
ity, x') for an a:: = 0.029 crystal at zero DC magnetic field 
and at pressures identified in the figure. Data taken with 
an excitation field of 1 Oe at a frequency of 1 Hz. b) Curie 
temperature, Tc, as defined by the temperature of the peak 
susceptibility in frame a, of our x = 0.029 crystal as a function 
of P. Also shown are the Weiss temperature, Qw, and Curie 
constant derived from fits of the Curie- Weiss form to the real 
part of the AC susceptibility shown in frame a. Note: Qw 
and the Curie constant were determined only for pressures 
where sufficient T > 15 K x' data (not shown) was available. 



C. The Specific Heat and Entropy 

The specific heat, C of our crystals was measured to 
characterize the changes to the thermodynamic proper- 
ties that occur upon Co doping FeS2. As can be seen 
in Fig. [HI where C is plotted for 7 representative crys- 
tals, there is little change in C{T) with Co doping be- 
tween 50 and 300 K. In frame b of the figure we con- 
centrate our attention on the low temperature specific 
heat where changes in the electronic contributions are 
apparent. Here we plot the data in the standard manner 
for characterizing the electronic contributions, C{T)/T 
vs. T^, since in standard paramagnetic metals the spe- 
cific heat is expected to be accurately described by a 
C{T)/T = 7-1- /3T^ form in the temperature range dis- 
played. In this equation 7 represents the coefficient of 
electronic specific heat and is proportional to the density 
of electronic states, while (3 parameterizes the phonon 



FIG. 9: (Color online) Specific heat, a) The temperature, 
T, dependence of the specific heat, C for seven representative 
Fei_a;Co3;S2 crystals with x identified in the figure, b) Specific 
heat divided by temperature, C/T, plotted as a function of 
for temperatures between 2 and 70 K. Symbols are the 
same as in frame a. 

We expand the low-T C{T)/T data in Fig. [TO] to high- 
light the changes that occur to the electronic contribu- 
tions with increasing x. In this figure the changes to 
7 with X above 1 K are made clearer along with the 
large increases in C{T)/T that occur below 1 K in the 
Co doped samples. The increases to 7 with x are very 
large considering that Hall effect measurements estimate 
the carrier concentrations are only 10 to 30% of the Co 
concentration in our samples^^. Such a large density of 
electronic states is commonly interpreted as a large elec- 
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tronic effective mass. Our estimates for the carrier effec- 
tive mass, m* , from the C{T)/T data above 5 K range 
from 4 times the bare electron mass, mg, to over 40 mg 
for samples with x just beyond Xc- If we were to con- 
sider the large upturn in C{T)/T below 1 K to be elec- 
tronic in origin, our estimates of m* increase to values 
between 120 and 250 mg at 120 mK. Heavy mass car- 
riers are rare in transition metal compounds, however, 
several Fe-based semiconducting, or semimetallic, mate- 
rials have been reported to have substantial carrier mass 
enhancemen t "^"i"^^ . Crystals with a; < Xg, as determined 
from the AC magnetic susceptibility, display a C{T)/T 
that increases continuously down to the lowest measure- 
ment temperatures (0.1 K). Those with x > Xc display a 
maximum in C (T) /T at a temperature somewhat below 
Tg as determined from the temperature of the maximum 
in x'- 



the low temperature C{T)/T while increasing slightly 
C{T)/T above 1 K. This behavior is qualitatively similar 
to the changes that occur upon application of a magnetic 
field in simple doped semiconductors such as Si:P'*^. The 
low-T specific heat of Si:P has been modeled in terms of 
local moments induced by the doping that are thought 
to form singlets such that the ground state of the system 
has zero net moment^i^. The effect of disorder is to 
create a broad distribution of interactions between the 
local moments. As a result, C{T)/T increases as a small 
power-law at low temperatures^. However, we do not 
believe such a description to be relevant in the case of 
Fei_j;Coj;S2 since our crystals have a nearly ferromag- 
netic, or ferromagnetically ordered, ground state. In- 
stead, we attribute the unusual temperature dependence 
of C{T)/T as being due to the interaction of a random 
distribution of local moments together with the interac- 
tion of the local moments and conduction electrons. As 
the magnetic field is increased beyond the level of the 
interactions, C{T)/T resembles a sum of contributions 
from a Fermi gas of electrons (a constant 7) and a sep- 
arate contribution from effectively noninteracting mag- 
netic moments. The contribution from non-interacting 
moments in a magnetic field is commonly referred to as 
a Schottky anomaly^ and is demonstrated in the figure 
by the dashed lines. We note that above 1 K the specific 
heat in zero field is well described by this Schottky form 
added to a temperature independent 7 and the phonon 
contribution taken as C(T) /T of the nominally pure FeS2 
crystal. 



FIG. 10: (Color online) Low temperature specific heat, a) 
The specific heat divided by temperature, C/T, plotted as 
a function of below 10 K for nine representative crystals 
with X identified in the figure, b) Specific heat divided by 
temperature, C/T, plotted as a function of T on a logarithmic 
scale. Symbols are the same as in frame a. Dashed lines are 
fits of the form aT"'^'^ to the data with an ac of -0.45±0.05 
for the X = 0.0007, -0.70±0.05 for the x = 0.005, -0.60±0.05 
for the X = 0.007, -0.70±0.05 for the x = 0.024, and -0.50± 
0.05 for the x = 0.045 crystals. Note that the temperature 
range for the fits was limited to T > IK for the x — 0.024 
and X = 0.045 crystals. 



In order to characterize the specific heat further, we 
have measured the magnetic field dependence of C{T)/T 
as shown in Fig. [TT] for 4 of our single crystals. For all 
four of the samples application of magnetic field decreases 
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FIG. 11: (Color online) Magnetic field dependence of the low 
temperature specific heat, a) The magnetic field, H, and 
temperature dependence of the specific heat divided by tem- 
perature, C/T for an a; = 0.024 crystal at fields identified in 
the figure. Dashed line is a fit of the Schottky forrr4^ for a 
collection of free magnetic moments in a magnetic field added 
to the phonon contribution to the specific heat determined by 
C of a nominally pure single crystal to the H = C/T data 
for T > 2 K. This fitting procedure was carried out with the 
quantity ng^J{J + 1) held constant at the value determined 
by fits of the Curie- Weiss form to the real part of the AC 
magnetic susceptibility measured on the same crystal. The 
fitting procedure resulted in best fit parameters of the effec- 
tive field, Heff ~ 1.4 T, a magnetic moment, J = 2.2, and a 
density of magnetic moments, = 9.5 x 10^^ cm"'', b) The 
H and T dependence of C/T of an a: = 0.007 single crystal 
at fields identified in the figure. Dashed line is the a fit of 
the same form as in frame a to the H = data. Best fit 
parameters were Heff = 1.8 T, J = 2.6, and A = 2.0 x 10^^ 
cm"^. c) The H and T dependence of C/T of x = 0.005 and 
X = 0.0007 single crystals as identified in the figure at fields 
identified in the figure. Dashed line is the a fit of the same 
form as in frame a to the x — 0.005, H = data. Best fit 
parameters were Heff = 1.88 T, J = 2.4, and A = 1.6 x 10^^ 
cm~^. 
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We have checked that our conclusion of formation of 
ferromagnetic cluster of magnetic moments at low-T is 
consistent with our C{T) measurements by quantita- 
tively comparing C{T)/T and xdc- As we stated above, 
the large Curie constants seen in Fig. [4] are likely caused 
by the fluctuations of clusters of magnetic moments and 
both the specific heat and the magnetic susceptibility are 
sensitive to the density and size of the fluctuating mo- 
ments. Therefore, we have fit the specific heat at T > 2 
K as shown in Fig.[TT]with the Schottky forrti^ discussed 
above for a set of non-interacting magnetic moments sub- 
jected to an effective molecular field. This contribution 
was added to a T-independent term representing a Fermi 
liquid of conducting electrons and the specific heat of the 
nominally pure FeS2 crystal as a model of the contribu- 
tion to C(T) from phonons. For this fitting procedure we 
have held the the quantity ng^J{J + 1) constant at the 
value determined by Xdc{T) and calculated the best fit 
values of the effective molecular field, H^ff, J, and n. In 
Fig.fT^ we present the results, plotting the best fit values 
of J, which varied from less than 3 to 6 across our range 
of X, and the parameter 2NJ/x where N = n/npu and 
npu is the density of Fei_2;Co2;S2 formula units. We plot 
the quantity 2NJ/x to check our assumption that each 
Co dopant donates a local magnetic moment of s = 1/2 
to FeS2. A value of 2NJ/x = 1 corresponds to agree- 
ment with this assumption and although there is signifi- 
cant scatter in Fig. [T2b. there is general agreement. The 
best fit values of i?e// varied between 1.4 and 2.6 T. We 
conclude that both x' and C{T)/T above 2 K are con- 
sistent with fluctuating magnetic moments whose size is 
several times larger than the magnetic moment of a single 
non-interacting Co dopant in an FeS2 lattice. 



C(T) prior to performing the integrations and the re- 
sult of this procedure is shown in Fig. [T31 Here S ap- 
proaches a;i?ln(2) only for temperatures well above 10 
K. Thus, there is significant entropy missing at 10 K 
consistent with the idea that clusters of spins form at 
temperatures of this order. Taking the entropy at 7 K 
as a benchmark allows us to estimate the average clus- 
ter size at this temperature. We make use of the form 
S = {x/k) R ln{2{k/2) + 1) where k is the average clus- 
ter size of spin 1/2 moments, presumably localized on 
Co sites, to make the estimate. We find average cluster 
sizes that range from about 3 to 5 for the samples where 
data below 1 K are available. If we assume an expo- 
nentially decaying probability for clusters containing TV 
spin-1/2 Co impurities, we estimate that about 10% of 
the cluster sizes will contain 10 or more spins at this T. 
The reduction of entropy below xR In 2 with cooling ev- 
ident in Fig. 1131 implies long ranged interaction between 
the local magnetic moments induced by Co substitution. 
For these Co concentration levels, and assuming that Co 
ions are substituted randomly for Fe atoms in the FeS2 
crystal structure, we calculate that the average distance 
between impurity sites to be at the level of third nearest 
neighbors. Thus we estimate the range of interactions be- 
tween the magnetic moments to be at least 1 nm in order 
to produce the reduction of entropy we observe below 10 
K, as well as the formation of magnetic clusters apparent 
in the magnetic susceptibility to somewhat higher T. 



FIG. 12: (Color online) Magnetic moments from specific heat 
and susceptibility analysis. Magnetic moment, J, and the 
quantity 2NJ/x where is the density of magnetic moments 
of size J per Fei_iCoiS2 formula unit determined from fits 
to the DC susceptibility and the specific heat above T = 2 K 
as described in text. 

The conclusion of magnetic cluster formation was fur- 
ther explored by calculating the entropy, S{T), from C 
according to the relation C{T')/T'dT' = S{T) for the 
small number of crystals that we have measured below 
1 K. We have subtracted the phonon contributions to 



FIG. 13: (Color online) Entropy determined from the specific 
heat, a) The temperature, T, dependence of the entropy, S, 
determined by numerical integration of the C/T data for five 
representative single crystals with stoichiometry's identified 
in the figure. 

In addition, the changes to the entropy that occur with 
magnetic field can be viewed in Fig.[T31 The sensitivity of 
C{T)/T and S to magnetic fields confirms the magnetic 
origin of the specific heat in this temperature range. For 
very small x, x = 7 x 10"'* for example, increases in S are 
observed below 2 K for iJ = 1 T, a sign that there is still 
a significant contribution to the entropy that we are not 



12 



accessing because it lies below our lowest measurement 
temperature. Thus, for small x, there is a significant 
density of magnetic moments that interact with energy 
scales of less than 100 mK. At larger x, x > 0.007, a 
suppression of S with field is observed from which we 
conclude that there are very few magnetic moments with 
such small energies of interaction. This is consistent with 
our conclusion that there is significant cluster formation 
in these samples at temperatures above those accessed in 
Fig. [T3] and for our x = 0.024 sample, a peak in x' at 
0.5 K that we associate with a disordered ferromagnetic 
transition. 



FIG. 14: (Color online) Field dependence of the entropy de- 
termined from the specific heat. The temperature, T, depen- 
dence of the entropy, S, determined by numerical integration 
of the C/T data at magnetic fields, H identified on the figure 
for a) an a; = 7 X 10"*, b) an a; = 0.007, and c) an a; = 0.024 
crystal. 



atically above the solid line representing the Curie con- 
stant expected when each Co dopant donates a fluctuat- 
ing J — 1/2 magnetic moment. If we focus our attention 
on those samples closest to Xc, we find Curie constants 
1.5 to 1.75 times larger than the simple estimate based 
on the expectation that each Co dopant adds a single 
spin-1/2 moment and the Curie law, x{T) = CC /T with 
CC — g^^^nJ{J + l)/3fcs. This is an indication that 
FM ordered clusters of magnetic moments are forming 
below 100 K, where the Curie behavior becomes evident 
for this dopant concentration. With the assumption that 
the magnetic moments are coupled ferromagnetically, the 
size of the Curie constant would indicate that on average 
each cluster contained between 2.5 and 3.2 spin-1/2 mag- 
netic moments for x x^. In addition, we have checked 
that this estimate is consistent with C{T) above 2 K 
where the data can be interpreted in terms of a Schottky 
anomaly associated with entropy of fluctuating magnetic 
moments. This estimate of the average size of fluctuat- 
ing magnetic clusters is further supported by estimates 
made from the entropy determined from the integral of 
C{T)/T (Fig. [131). 

Motivated by theoretical predictions about the impor- 
tance of rare regions in disordered magnetic systema^i^S, 
we have carefully examined the temperature and field de- 
pendent behavior of the magnetization and specific heat 
of our samples with x in proximity to Xc to search for 
such effects. We are not disappointed in that unusual 
behavior is discovered. The most clear indication of the 
importance of rare large clusters of magnetic moments 
comes from the temperature dependence of the specific 
heat below 1 K shown in Fig. [TOb. Here C{T)/T is seen 
to increase without any suggestion of saturation down 
to 0.1 K for our samples closest to Xc, x = 0.005 and 
0.007. As indicated by the quality of the fits of a power- 
law form, represented by the dashed lines in the figure, 
as well as the linearity of the data when plotted on a 
log-log scale in Fig. [12] a, these data are well represented 
by a C(r) oc T"^"! form with ac between 0.3 and 0.45 
for over a decade in T. Similar power-law behavior is 
observed in x' over the same temperature range for the 
same crystals as can be seen in Fig. [T5]b. Here x' is well 
fit by a T"^~^ form between 0.8 and 10 K with be- 
tween 0.2 and 0.25 for these samples. We note, however, 
that x', unlike C{T)/T tends to saturate below 0.5 K. 

We make a direct comparison of the real part of the 
AC magnetic susceptibility with C{T)/T in Fig. [TBI where 
the Wilson ratio, defined as 



IV. GRIFFITHS BEHAVIOR NEAR THE 
CRITICAL CONCENTRATION 

As we noted in our previous discussion of both the 
AC and DC magnetic susceptibility, the Curie constants 
determined from simple fits to the temperature depen- 
dence were larger than is to be expected for paramag- 
netically fluctuating spin-1/2 magnetic moments. That 
is, the Curie constants displayed in Fig. [3| lie system- 



kl x7 C{T)/T 



(1) 



is shown. In this figure we observe that the two sam- 
ples with X close to Xc have a Wilson ratio somewhat 
larger than 10. This value is similar to that observed near 
the critical point for magnetic ordering in CePdi-^jRha, 
which is thought to enter a novel Kondo-cluster-glass 
phase and has properties compatible with the quantum 
Griffiths phase scenario^. Estimates for the average sizes 
of the fluctuating magnetic moments from the Wilson ra- 
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FIG. 16: (Color online) Wilson Ratio. The Wilson ratio, 
■K^k^x I'^AC /T , for crystals with x identified in the figure. 



FIG. 15: (Color online) Low temperature specific heat and 
magnetic susceptibility, a) The temperature dependence of 
the specific heat divided by temperature, C /T for several rep- 
resentative crystals with x identified in the figure. Dashed 
lines are fits to C/T to a power-law form, C{T)/T oc T°'c-^ 
as in Fig. 1101 b) real part of the AC magnetic susceptibility, 
x' for the same crystals as in frame a. Symbols the same as 
in frame a. dashed dotted lines are fit of a power-law form 
to x', X' « T"'^'^ with exponents, of 0.25 ± 0.1 for the 
X = 0.005, 0.22 ± 0.1 for the x = 0.007 and -0.3 ± 0.2 for the 
X = 0.024 crystals. 



scatter in the values is apparent along with the tendency 
for larger a values to be found in specific heat measure- 
ments than in magnetic susceptibility measurements. 



tio of Fei_2: 00^:82 yield about 3 to 4 spins consistent with 
our other estimatesiS. We note that the Wilson ratio at 
the lowest temperatures is not constant, but instead re- 
tains a slight temperature dependence reflecting the ob- 
servation that at low-T C{T)/T continues to display a 
power-law temperature dependence whereas x' tends to 
saturate. We note that predictions of the Wilson Ratio 
in Griffiths phases models indicate a weak, logarithmic, 
temperature dependences^. For our samples with x > Xc 
the Wilson ratio becomes very large exceeding 1000 near 
Tc indicating a strong ferromagnetic tendency. 

In addition, our finite-i/ DC magnetic susceptibility 
measurements indicate a small power-law temperature 
dependence below 10 K as can be seen in Fig. [iTl Al- 
though the temperature range is restricted, for samples 
with a; < a fit of the form xdc oc T"""^"^ results in 
small, between 0.4 and 0.05, values of aoc- We note that 
these data were taken in significant DC magnetic fields of 
between 1 kG and 1 T in order to produce reasonable sig- 
nal sizes for our SQUID magnetometery measurements. 
Even so, we find similar small power-law behavior with 
the tendency for aoc to fall toward zero near Xc- A 
compilation of the power-laws determined from these dif- 
ferent measurements is presented in Fig. [TH] where some 



FIG. 17: (Color online) Power-law temperature dependence 
of the DC magnetic susceptibility. The DC magnetic suscep- 
tibility, Xdc, vs. temperature, T for crystals with x < Xc 
as identified in the figure. Data is the same as that shown in 
Fig.[l Lines are fits to the form xdc = XPauU + CC /T"'"^-^ 
as described in the text. 

Our data and analysis outlined above have led us to the 
conclusion that ferromagnetic clusters form at T > Tc- 
The sensitivity of our susceptibility measurements to 
small increases in pressure is one of the consequences of 
these weakly interacting ferromagnetically coupled clus- 
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FIG. 18: (Color online) Results of power-law fits to the mag- 
netic susceptibility and specific heat, (a) The best-fit expo- 
nents from fits of simple power-law forms to the AC and DC 
susceptibility and specific heat, C, temperature, T, and mag- 
netic field, H, dependence above the Curie temperature of 
Fei_^Co^S2 for 3 x 10"" < a; < 0.055. Here aoc is the 
temperature exponent from fits of the form Xdc = cT°'°'^~^ 
to the DC magnetic susceptibility data below 10 K, ^ deter- 
mined from a scaling of M{T, H) data described in the text, 
A taken from fits of Eq. [2] to T-dependence of the real-part 
of the AC susceptibility, x', from fits of x' to the form 

X' = cT"^-^ between 0.8 and 10 K, ac taken from fits of 
the form C/T = cT"^''^ to the T-dependence of the spe- 
cific heat, Oh taken from fits of the form C'/T — cH°'"~^ to 
the //-dependence of C/T, and qa/ is taken from fits of the 
form x' = cH°'-^ to the magnetic field, H, dependence of x' ■ 
Dashed line is a guide-to-the-eye. (b) Plot of the temperature 
characterizing the maximum interaction between magnetic 
moments, taken from fits of Eq. [2] to the T-dependence 
of x' for T > n. 



T-dependence suggests a magnetization that is caused by 
independent fluctuators having a power-law local energy 
distribution, P(A) ^ A^^^, where A is the local energy 
scale of the fluctuator. The form of the magnetization for 
freely fluctuating paramagnetic moments is well known to 
follow the Curie law M(T, H) = ng^iBJBj{g^iBH/kBT), 
where Bj is the Brillouin function'*^ and thus scales as 
M{T, H) oc F{H/T). In the case of the power-law energy 
distribution the magnetization has been shown to scale 
as M(T, H) = H^F{H/T)^. As can be seen in Fig. [Mi, 
which displays M{T,H) between 1.8 and 10 K for our 
two crystals with x closest to Xc, the the standard scal- 
ing form does not scale our data satisfactorily. This is to 
be expected for a system on the verge of magnetic order- 
ing where magnetic moments have significant interaction. 
However, by simply allowing for a small value of ^ we can 
scale all of our H and T dependent data reasonably well 
as can be seen in frame b of the figure = 0.15). We 
note that ^ is smaller than the values we found for ac and 
a DC 7 but still distinct from zero. We note that this form 
for M does not hold well for x > 0.04 since M{H) in- 
creases sharply at low fields thus resembling more closely 
M{H) of ferromagnetic materials. 



ters. In addition, we observe small power-law dependen- 
cies of C{T)/T, x\ and xdc at low-T with changes in 
temperature as was demonstrated in Fig. [15] for in 
Figs. [TOl and [T5l for the specific heat, and is apparent in 
Fig. [T7| for Xdc- In each case these quantities are ob- 
served to vary as T"^^ with < a < 0.6. All of these 
power-laws have similar values although there is signifi- 
cant scatter apparent. However, the similarities between 
these power- laws, aM, oldc, and ac for individual sam- 
ples implies that a single physical mechanism may be 
responsible. 

A clue as to the identity of the physical mechanism 
can be found in the dependence of the magnetization on 
temperature and magnetic field, M{T,H). As we have 
noted above, the magnetic susceptibility does not con- 
form to Curie- Weiss dependence, instead being charac- 
terized by a smaller power-law divergence at low-T. This 



FIG. 19: (Color online) Scaling plot of the Magnetization, a) 
Plot of the magnetization, M, of Fei_a;Coa;S2 crystals with 
X — 0.007 (left side scale) and x = 0.012 (right side scale) as 
a function of magnetic field, H, divided by temperature T, 
(H /T) to display the quality of the standard paramagnetic 
scaling form M{T,H) oc F{H/T). Constant temperature, T, 
sweeps are shown for T's indicated in the frame and in frame b 
for 1.8 < T < 10 K for < 5 T. b) Plot of the same A/(T, H) 
data as in frame a, divided by as a function of H/T. 
Reasonable scaling of the data is achieved for ^ = 0.15 ± 0.05. 

In addition to this unusual field and temperature de- 
pendence of M, our specific heat measurements also re- 
veal a strong sensitivity to magnetic fields as we demon- 
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strated in Fig. [TT] Although we found that above 2 K 
the T and H dependence of C/T appears consistent with 
a Schottky anomaly form, C/T below 1 K is not, even 
when a 3 T magnetic field is applied. This is emphasized 
by the magnetic field dependence of the low-T specific 
heat which is explicitly displayed for T = 120 mK Fig.[20l 
Here moderately sized fields are seen to suppress C (T) /T 
dramatically. The power-law like decrease in C{T)/T is 
displayed in the inset where the lines are power-laws in 
field, C(T)/T oc H""-'^ with = for the x = 0.005 
crystal and 0.4 for the x = 0.024 crystal. We note that 
for a Schottky anomaly C{T)/T decreases exponentially 
in H for fields where gfJ-sH > ksT which is inconsistent 
with our datc4^. 



FIG. 20: (Color online) Magnetic field dependence of the spe- 
cific heat at T = 0.12 K. The magnetic field, H, dependence 
of the specific heat divided by temperature, C/T for several 
crystals with stoichiometry's identified in the figure. Inset: 
The same C{T, H) /T data plotted with logarithmic axes. The 
lines are fits of a power-law form to C /T data for the x — 0.005 
and X — 0.024 crystals. Best fit values of the exponents are 
-1 for the X = 0.005 crystal and -0.6 for the x = 0.024 
crystal. 

Although there are many physical realizations of this 
scaling form for the magnetization we observe (indepen- 
dent fluctuators with a power-law energy distribution), 
and small power-law T dependent magnetic susceptibili- 
ties have been known to describe doped semiconductors 
very near insulator-to-metal transitions, our identifica- 
tion of magnetic clusters forming in this doped semicon- 
ducting system presents an obvious candidate explana- 
tion for all of the data presented thus far. We posit 
that our data lead naturally to the conclusion that Grif- 
fiths phase formation is the most likely physical mecha- 
nism responsible for the unusual properties we measure 
for samples in proximity to the critical concentration for 
magnetic ordering. 

Magnetic Griffiths phenomena have been used to de- 
scribe disordered magnetic systems theoretically in the 
case where the disorder is sufficient to cause inhomogene- 
ity in the formation of the magnetic ground state^^. In 



the vicinity of the quantum critical point, droplets of the 
magnetically ordered phase are thought to form due to 
the inhomogeneities inherent to a chemically substituted 
material. These droplets reside within the paramagnetic 
phase and can be exponentially rare if the statistics are 
Poissonian so that P{V) ^ e~'^^ , where P is the probabil- 
ity of having a droplet of volume V and c is related to the 
strength of the disorder^^. Since Fci-^: 00^:82 has a ferro- 
magnetic ground state for x > Xc, the droplets most likely 
carry a net magnetic moment simply proportional to the 
size of the droplet. These magnetic moments will have a 
finite probability of tunneling to nearby excited magne- 
tization states. The tunneling rate for such a droplet is 
predicted to be exponential in the number of spins that 
form the droplet, such that A ~ ujQe~^^ where b is a 
constant related to the microscopic tunneling mechanism, 
and ujQ is a frequency cutofS^. Since the tunneling rate is 
proportional to the energy splitting of the different mag- 
netization states the distribution of energy splittings can 
be written as P(A) ~ / dVP{V)S[A - cj^e-''^] - A^-\ 
where ^ = c/b, giving us our power-law distribution of en- 
ergy levels suggested by the scaling of our M{T, H) data. 
With this form for the energy splittings, power-law forms 
for the magnetic susceptibility and specific heat follow 
by considering all clusters with A > fc^T to be frozen 
such that they do not contribute to the susceptibility or 
entropyiS,. Thus, a modified Curie form x'{T) oc n{T)/T, 
where 'n{T) is the density of clusters with A < fc^T is 
predicte d^^i^° . n(T) can be easily calculated from P(A) 

as n{T) oc P{A)dA oc and thus, x'{T) oc T«-i. 
In addition, ^(r) oc n(r)i?ln (2J -1- 1) oc T" and since 
C{T)/T = dS/dT we have C{T)/T oc T^"". There- 
fore, the Griffiths phase phenomenology leads to a form 
of the magnetization, susceptibility, and, C{T)/T similar 
to that found in our experiments. 

There has been some discussion in the literature about 
the ability of magnetic clusters to tunnel into different 
magnetization states in the presence of itinerant elec- 
trons. The presence of the conduction electrons which 
couple to the local magnetic moments via the Kondo ef- 
fect cause dissipative effects on the tunneling of magnetic 
clusters. The dissipation limits the size of clusters al- 
lowed to tunnel into opposing magnetization states thus 
cutting off the singular behavior in the models^^. Mil- 
lis et al. argue that in heavy fermion metals the the 
electrons are sufficiently well coupled to local magnetic 
moments, with Kondo temperatures of order 100 K, so 
that Griffiths phases are unlikely to be observed in these 
systems^S. The idea being that the dissipation is relevant 
at the scale of the Kondo temperature, thus the temper- 
ature region where Griffiths phases may be observed is 
very small. However, Castro Neto and Jones argue that 
in the region of magnetic clusters the coupling of conduc- 
tion electrons to the magnetic moments is not important 
since in these regions the electron density of states is not 
renormalized by the Kondo effec t^^i^^ . 

In Fei_a;CoxS2 the relevant energy scales are such that 
Griffiths phases may be observed over a much wider tem- 
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perature scale. Here, the Kondo coupling scale is only of 
order 2 K in the region of the critical Co concentration 
for ferromagnetism. The Kondo scale is likely so small 
because the density of electron states, in the approxima- 
tion of a spherical electron pocket, is proportional to n'J'^ 
where rie is the electron density which is more than 2 or- 
ders of magnitude smaller, for the x range investigated 
here, than in most heavy fermion metals^^. In addition, 
the disorder inherent in a doped semiconducting system 
results in a high scattering rate of the electrons further 
reducing their ability to screen local moments. There- 
fore, we expect the dissipation due to a conducting gas 
of electrons to be a much smaller effect than in the heavy 
fermion metals where electron densities are much larger 
and the strong Kondo coupling renormalizes the electron 
DOS to very high levels. In addition, the relevant en- 
ergy scale for local moment interactions is commonly es- 
timated by the Curie temperature of the clean magnetic 
system without dilution. In this case the clean ferromag- 
net is C0S2 and this sets the largest energy scale for local 
moment interactions at 120 K. Estimates for the average 
interaction energy of local moments include the Weiss 
temperature which we have measured via the magnetic 
susceptibility and is displayed in Fig.[4l From the average 
and maximum local moment interaction scales compared 
with the Kondo temperature scale in our materials, we 
conclude that dissipative effects of the conduction elec- 
tron gas should be much smaller in Yqx^^Cox'S)'! than in 
heavy fermion metals. However, rounding of the tran- 
sition due to Kondo coupling of the local moments and 
the conducting electrons may be relevant at the lowest 
temperatures. 



displayed for frequencies between 1 and 1000 Hz. Here, 
we observe significant decreases in the magnitude of both 
x' and x", very small changes to the temperature of the 
maximum in x' (at most 0.1 K increase), and a larger 
increase of the temperature of the maximum in x" with 
measurement frequency. We observe no indications of the 
onset of significant skin effects and eddy current absorp- 
tion which cause a steady increase in the out-of-phase 
susceptibility (x") with / despite the large out-of-phase 
component of the signal. 

The frequency dependence we measure is similar to the 
frequency dependent changes that occur in the archety- 
pal spin glass CuMn^ with a relative shift of per 
decade of measurement frequency of less than 1%. How- 
ever there are with two important differences with typical 
spin glass behavior that draw our attention. First, we ob- 
serve much larger changes in the magnitudes of both x' 
and x" with frequency than is typical for a metallic spin 
glass. Second, the temperature range over which we ob- 
serve a suppression of x' is from at least 0.5 T^. to ~ 2Tc, 
whereas for CuMn of similar dopant concentrations there 
is no measurable /-dependence of x' for reduced temper- 
atures, (T — Tf)/Tf > 0.05, where Tf is the spin glass 
freezing temperature^. 



V. FINITE TEMPERATURE MAGNETIC 
TRANSITIONS; x > 

The figures and discussion above give a general descrip- 
tion of Fei_a;Co3;S2 as it evolves from a diamagnetic in- 
sulator into a ferromagnetic metal. These data establish 
a magnetic transition that can be observed via the AC 
magnetic susceptibility down to very low temperatures at 
very small x . It is clear from our data that magnetic clus- 
ters form above Tc and an inhomogeneous magnetic state 
below Tc- In this section we explore in more detail the 
behavior of our Fei_2;Co2;S2 crystals having a magnetic 
phase transition, or perhaps a magnetic glass freezing, 
at finite temperatures. We examine the frequency and 
magnetic field dependence of the magnetic susceptibil- 
ity both above and below the ordering to compare with 
typical spin glasses and disordered FMs to gain insight 
on how such a low density disordered metal establishes a 
magnetic ground state. 

To this end we have carefully measured the /- 
dependence of x' and x" of a few of our samples with 
X > Xcto probe the low frequency dynamics of the highly 
paramagnetic state above Tc as well as the magnetic state 
below Tc- In Fig.[2T]x' and x" of our x = 0.045 sample is 



FIG. 21: (Color online) Frequency dependence of the AC sus- 
ceptibility. The frequency and temperature dependence of 
the (a) real part, x'j and (b) imaginary part, x", of the AC 
magnetic susceptibility at zero DC field for several x = 0.045 
single crystals. Measurement frequencies were varied between 
1 and 1000 Hz as identified in the figure. The excitation field 
was 1 Oe. 

The frequency dependence of x' and x" for several 
temperatures is explicitly demonstrated in Fig. [21] for 
temperatures near Tc- While x' decreases monotonically 
with frequency x" is characterized by a broad peak which 
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moves through the frequency window of our measurement 
from low / to high / with temperature. At T < T^, x'if) 
is well described by a logarithmic dependence, as demon- 
strated in the figure by the dashed lines, with an extrapo- 
lated zero crossing between lO'* and 10"'^'* Hz. This depen- 
dence has also been observed in insulating spin glasses 
such as Ising-like LiHo^^Yi-j^F^^. In addition, x" ti^is 
very little /-dependence in the same temperature range 
so that the dynamic response of the system below Tc is 
much like that of a spin glass where fluctuations occur 
on all long time scales?-. Thus, the magnetic dynamics 
at these low frequencies resemble that of common spin 
glass systems, with the added feature that the time de- 
pendence extends to temperatures far above Tc- 



a shoulder 7K as well. In addition x" has two max- 
ima, at 4.5 K and at ~ 7 K. The differences we observe 
in Fig. [221 with application of very small fields, as well 
as the appearance of structure in the T-sweeps when the 
field is carefully zeroed are extraordinary. We interpret 
this sensitivity as an indication of enormous degeneracy, 
or near degeneracy, of magnetic moment configurations 
in these crystals. We observe that DC fields as small 
as 10 Oe have a significant effect, decreasing the peak 
value of x' by a factor of 3.5 and x" by more than a 
factor of 30. The reason for the marked difference be- 
tween the AC and DC magnetic susceptibilities that we 
noted above is now quite clear; the DC magnetic suscep- 
tibility was measured at fields between 50 and Iff* Oe 
applied (DC) field, whereas x' shows dramatic changes 
between and 200 Oe in the range 1.8 < T < 20K for 
all samples measured with AC susceptibility techniques 
at finite H. Therefore, the AC susceptibility data repre- 
sent a clearer impression of the H — magnetic state of 
our Fei_a;Coa;S2 crystals. In addition, the sensitivity to 
frequency and magnetic fields indicates a complex energy 
landscape with nearly degenerate energy minima for the 
magnetic moments. Our main conclusion from the broad 
temperature and field dependence of M and C in sections 
nil A[ HITBI and llll CI that clusters of magnetic moments 
form at r >> Tc, is consistent with the conclusion of a 
spin-glass-like sensitivity to measurement / and applied 
magnetic fields over a broad T-range about Tc- 



FIG. 22: (Color online) Frequency dependence of the AC sus- 
ceptibility. The frequency dependence of the real part, x', (a) 
and imaginary part, x", (b) of the AC magnetic susceptibility 
at zero DC field for several x = 0.045 single crystals. Temper- 
atures shown are 1.8 K (red diamonds), 3 K (orange bullets), 
4 K (green filled squares), 6 K (blue boxes), and 9 K (violet 
triangles). The excitation field was 1 Oe. Dashed black lines 
in frame a are fits to a logarithmic frequency dependence of 

x'. 

In addition to the enhanced sensitivity to measurement 
frequency near Tc our crystals display an enhanced sensi- 
tivity to external magnetic fields in the same temperature 
range. This is demonstrated in Fig. [521 for an a; = 0.045 
sample consisting of several small crystals. We note that 
the data displayed in Fig. [52b was previously publishediS, 
and is reproduced here for completeness and compari- 
son of x' s-nd x" ■ In this experiment we have carefully 
minimized the DC magnetic field for our H = scan 
which reveals not only the x' maximum near 4.5 K but 



FIG. 23: (Color online) Dramatic field dependence of the AC 
magnetic susceptibility. The magnetic field, H, and tempera- 
ture, T dependence of (a) the real part, x'l and (b) imaginary 
part, x"i of the AC susceptibility, for the same sample as in 
Fig. [21] at fields identified in the figure. Data taken with an 
excitation field of 1 Oe at a frequency of 1 Hz. 
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To aid in understanding the changes that occur to the 
magnetic response of this system just above the order- 
ing temperature, we have compared 1/x' taken at small 
DC magnetic fields as a function of T to the Curie- Weiss 
form as shown in Fig. [24] for x — 0.023 and x — 0.045. 
What is interesting here is that 1/x' falls below the 
Curie- Weiss behavior established at higher temperatures. 
That is, x' diverges more strongly than the Curie- Weiss 
behavior for temperatures substantially above T^. We 
note that 1/x' evolves toward the Curie- Weiss form, 
1/x' = [T - Q^)/CC, with the apphcation of small 
fields of order 100 Oe for x = 0.023 and only 10 Oe 
for X — 0.045. The strong increase of x' above Tc and 
is evidence for short ranged ferromagnetic order, or 
as we have noted above, ferromagnetically aligned clus- 
ters of magnetic moment formation^ We conclude from 
this form of x' that the ground state of this system is 
likely disordered ferromagnetic rather than a spin-glass 
state. In the bottom half of Fig. [M] we demonstrate that 
a power-law temperature dependent for m^^i^'^ , 

i/x'cx(r/rj-i)i-^ (2) 

with A ~ 1/2, describes x' for 1-5 to 2 orders of mag- 
nitude in reduced temperature above TJ". Here TJ" is an 
indication of the temperature scale where ferromagnetic 
clusters begin to form and is an indication of the max- 
imum interaction strength between magnetic moments. 
In contrast, Qyp is an average, or mean field interaction 
between magnetic moments. A compilation of A and TJ" 
values from fits of Eq. [2] to 1 /x' for a small number of 
crystals with x varying between 0.012 and 0.045 is in- 
cluded in Fig. [181 While TJ" tends to increase with x, A 
takes on values between 0.2 and 0.4 for x > Xc- 

Beyond the enhanced frequency and magnetic field de- 
pendence seen in Figs.[2I]and[231 there are features in our 
data that indicate that the magnetic state of Fei_a;Coa;S2 
is far more interesting than that of a metallic spin glass, 
particularly for temperatures above Tc- For example, in 
Fig. [25|we plot the non- linear AC susceptibility, x'{H) 
and x"{H) at temperatures near Tc showing explicitly 
the changes that occur with H at constant T. Although 
this figure shows a broad, continuous, decrease of both 
x' and x" with H for T < Tc, it also demonstrates a very 
sharp low field behavior for T > Tc- A mean-field anal- 
ysis of the non-linear susceptibility involves the assump- 
tion of analyticity around H — such that a power series 
expansion in even powers, ~x'{H) — xi+X3^^+X5-f'^^ + --- 
can be employed to describe x'{H)- The extraordinary 
field dependence of x' seen in Fig.[25|cannot be described 
accurately by this form even over a reduced field range 
of only 10 Oe for T > Tc. Thus the randomness due 
to the chemical substitutional disorder has modified the 
field response of this system above Tc such that it is no 
longer analytic near H = 0. 

This form of x'{H) is similar to that found in (body- 
centered tetragonal) LiHoa;Yi_a;F4 where the magnetic 
moments associated with the Ho ions strongly align along 
the c-axis and couple to each other via the dipolar 



FIG. 24: (Color online) Temperature dependence of the in- 
verse AC magnetic susceptibility, (a) The temperature, T, 
dependence of the inverse of the real part of the AC magnetic 
susceptibility, 1/x' for x = 0.023, at DC magnetic fields, H, 
identified in the figure. The black line is the best fit of the 
Curie- Weiss form to the H = data between 20 and 50 K 
with OyfF = 3.0 K and a Curie constant corresponding to a 
population of 1.5 x 10^^ cm""^ spin 1/2 or 0.061 spin 1/2 per 
Feo.977Coo.o23S2 formula unit. The blue line is a fit of Eq. [2] 
in the text to the data with a best fit value of A = 0.44 ± 0.05 
and Tc of 10.5 ± 0.5 K identified in the figure by the red 
arrow, (b) 1/x' of our x = 0.023 sample as a function of re- 
duced temperature, T/T^ — 1, with determined from the 
best fit of Eq. [2]to the data above 11.5 K. The Blue line is the 
same as in frame a of the figure, (c) T-dependence of 1/x' for 
X = 0.045, at i/'s identified in the figure. The black line is the 
best fit of the Curie- Weiss form to the H = data between 
25 and 100 K with Qp^ = 12.8 ± 0.5 K and a Curie constant 
corresponding to a population of 1.8 x lO'^^ cm^"^ spin 1/2 or 
0.072 spin 1/2 per Feo.955Coo.o45S2 formula unit. The blue 
line is a fit of Eq. [2] in the text to the data with a best fit 
value of A = 0.55 ± 0.05 and TJ of 24 ± 0.5 K identified in the 
figure by the red arrow, (d) 1/x' of our x = 0.045 sample as a 
function of T /T// — 1 with determined from the best fit of 
Eq.[2]to the data above 26.5 K. The Blue line is the same as 
in frame c of the figure. Data taken with an excitation field 
of 1 Oe at a frequency of 1 Hz. 



interaction^^'^'^. The ground state is ferromagnetic for 
X > 0.2 and a spin-glass for 0.1 < a; < 0. 2^^'^'^ . A mag- 
netic field applied transverse to the moments, Ht , causes 
a mixing of excited states with the ground state doublet. 
This suppresses the ordering such that a quantum crit- 
ical point [T = 0) has been accessed for a large range 
of X. For X < 0.5 Refs. IHilli find that x'{Ht 0) 
and x'{Hi ~» 0), where Hi is along the c-axis, both have 
a finite slope, dx' /dH, and thus cannot be expressed in 
a perturbative expansion in H . The interpretation was 
that the disorder induced field distribution contains rare 
large amplitude random fields that prevent a description 
of x{Ht) in a mean field way as — > 0. The authors con- 
clude that the non-analyticity at TJt = and T > Tc is a 
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between and 1 describes the data over at least a decade 
and 1/2 in field for T > T^, while for T <Tc& much more 
gentle field dependence is seen. The best fit values 
have been included in Fig. (TH] for two samples where we 
have sufficient data for an accurate determination of this 
parameter. In addition, x" can also be described in a 
power-law form with an exponent between -1.3 and -2.0. 
For fields between 50 and 100 Oe x" is seen to undergo a 
much steeper decrease so that above these fields it is con- 
sistent with zero being no larger than the backgrounds in 
our measurement. The small power-law form for x'(-ff) 
at r > Tc is consistent with the Griffiths phase hypoth- 
esis made on the basis of the temperature dependence 
above. 



FIG. 25: (Color online) Magnetic field dependence of the AC 
magnetic susceptibility. The magnetic field, H, and temper- 
ature dependence of (a) the real part, x', and (b) imaginary 
part, x", of the AC susceptibility, for the same single crystal 
as in Figs. I2lland l23l at temperatures identified in the figure. 
Data tafcen with an excitation field of 1 Oe at a frequency of 
1 Hz. 



manifestation of Griffith singularities^^. Although there 
are similarities between our system, Fei_2;Co2;S2, and 
LiHo2;Yi_a;F4, both are disordered ferromagnets that 
are derived by chemical substitution from an insulat- 
ing parent compound, there are significant differences 
which are important to point out. The most obvious 
is that Fei_a;Coa;S2 is a metal for x > 3 x 10~'', while 
LiHoj;Yi_a;F4 remains insulating for all x. In addition, 
LiIl0j.Yi_2;F4 is strongly Ising-like with the magnetic 
moments preferring to align along the c-axis. In con- 
trast, Fei_a;Coa;S2 is cubic and for x = 1 has been shown 
to have a very small cubic anisotropy of its magnetization 
so that it been described as Heisenberg-like^^. However, 
we find that both of these disordered magnets contain 
rare large amplitude fluctuations that are the hallmarks 
of Griffiths phase formation for T > Tc- Thus, not only 
do we find evidence for such phases in samples on the 
verge of magnetic ordering at T = 0, but also discover 
evidence for similar rare region effects at larger x, x > Xc 
above the ordering temperature. 

In contrast to mean-field forms, the magnetic field de- 
pendence of x' shown in Fig. [25] can instead be well de- 
scribed as a simple power-law in H at T > Tc as demon- 
strated in Figs. and [571 These figures demonstrate 
that a small power-law form, x' — XoH""~^ , with aM 



FIG. 26: (Color online) Power law magnetic field dependence 
of the AC susceptibility. The magnetic field, H dependence 
of (a) the real part, x', and (b) the imaginary part, x", of the 
AC magnetic susceptibility with logarithmic axis for the same 
X — 0.045 sample as in Figs. 1211 1231 and at temperatures 
identified in the figure. The lines are power-law fits to the 
data with best fit exponents of -0.6 (-2.1) at 4 K, -0.8 (-1.5) 
at 6 K, -0.9 (-1.7) at 8 K and -1.0 (-1.9) at 10 K for x' (x")- 
Data taken with an excitation field of 1 Oe at a frequency of 
1 Hz. 

Fig. [771 displays x' and x" at 1.8 K for 5 different sam- 
ples with X > Xc and demonstrates several features that 
are common. First we note that for these Co concen- 
trations Tc > 1.8 K and, as discussed above, x' is not 
strongly field dependent below about 100 Oe at this tem- 
perature. This is in contrast to the data shown above Tc, 
a small subset of which is shown in the figure for compar- 
ison. For larger fields, x' is seen to decrease as ff-0-9±o.i 
while x" is suppressed to within the background level of 
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zero above this same field scale. It is interesting to note 
that x' above 10"^ Oe is less dependent on x suggesting 
that here the magnetic energy level distribution is much 
more independent of x. 



phase phenomenologj^ii^i22iiS^ where the non-universal 
power-law dependent properties are the result of quan- 
tum mechanical tunneling of these rare regions to states 
with different magnetizations. 



VI. DISCUSSION AND CONCLUSIONS 



FIG. 27: (Color online) Magnetic field dependence of the AC 
susceptibility for various Co densities larger than critical con- 
centration. The magnetic field, H dependence of (a) the real 
part, x' 1 ^nd (b) the imaginary part, x", of the AC magnetic 
susceptibility at 1.8 K with logarithmic axis for Fei_i:Co2:S2 
for x's identified in the figure. Data at one T > Tc are shown 
for the X = 0.024 (5 K) and x = 0.045 (8 K) crystals to 
demonstrate the changes that occur to the low field behavior. 
The solid lines are power-law fits to the data with best fit 
power-law of -0.9 for the high field data at T = 1.8 K as well 
as for the data for x — 0.045 at 8 K. Dashed line is a 
dependence for comparison. Data taken with an excitation 
field of 1 Oe at a frequency of 1 Hz. 

Our presentation in this section has focused on the AC 
magnetic susceptibility of our Fei_a;Coa;S2 crystals with 
X > Xc which demonstrates that we have discovered sev- 
eral unusual aspects. This includes a highly frequency 
dependent x' ^^'^ x" which encompasses a temperature 
range both above and below Tc that is much larger than 
in prototypical spin-glass systems. We have found an ex- 
traordinary non-analytic field dependence to the AC sus- 
ceptibility near H = coupled with a small power law 
form to the temperature dependence. These aspects are 
indicative of an inhomogeneous magnetic system where 
rare ordered regions form at temperatures above the long- 
range ordering, or perhaps a glass-like freezing, temper- 
ature of the disordered system. These local ferromagnet- 
ically ordered regions are well described by the Griffiths 



In this paper we have described a set of experiments 
exploring the nucleation of a magnetic ground state in a 
carrier-doped nonmagnetic insulator. For this investiga- 
tion we chose the relatively simple diamagnetic insulator 
FeS2 as it allows for Co substitution for Fe without the 
formation of second phase o^'^d^d^dSiiT addition, end 
member of the series C0S2 is an itinerant fcrromagnet 
with a Curie temperature of 120 K-'^'*. We found that the 
Co substitution yields a small number of itinerant charge 
carriers^^ and magnetic moments that are likely localized 
to the Co impurity sites. A weakly coupled magnetic 
ground state, spin glass or disordered ferromagnetic, de- 
veloped for X > 0.007 ± 0.002 as evidenced by a peak 
in the temperature dependence of the magnetic suscep- 
tibility. As X is increased beyond Xc the magnetic state 
becomes more robust to temperature and magnetic field 
so that the FM state is more evident at fields and tem- 
peratures where common tests can easily be performed. 

In the paramagnetic regions of temperature and Co 
concentration, we found fluctuating magnetic moments 
whose size was significantly larger than that expected for 
spin-1/2 moments of individual Co dopants in an FeS2 
background. This led us to the conclusion that clusters 
of magnetic moments were forming upon cooling below 10 
K for samples with x < Xc and at T > Tc for sample with 
X > Xc- The formation of clusters of average size of 3 to 5 
J = 1/2 magnetic moments was also found to be consis- 
tent with the measured C (T) /T and the entropy derived 
from it. Because of the disorder and assumed random 
placement of Co impurities throughout our samples, it is 
likely that a distribution of cluster sizes results. This dis- 
tribution is likely to have a very long tail for large cluster 
sizes so that there may be rare regions of local magnetic 
order. Such a description fits in well with the sensitivity 
to magnetic fields that we observed in the susceptibility 
and magnetization of our crystals. 

We posit that it is these rare regions of incipient order 
that cause the unusual temperature and magnetic field 
dependencies that we measured in the magnetic suscep- 
tibility and specific heat. The consistent appearance of 
small power-laws in our analysis of both the low temper- 
ature behavior of samples in proximity to Xc and at tem- 
peratures just above Tc for samples with x > Xc lead us 
to conclude that Griffiths phase physics is the most likely 
mechanism determining the static and dynamic proper- 
ties of this systeroi^iSS. In this scenario large rare re- 
gions of order appear in disordered magnets at tempera- 
tures below the critical temperature of the system with- 
out disorder. These ordered regions which act as single 
large magnetic moments are subject to tunneling events 
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where the magnetic moment reverses. When a distri- 
bution of cluster sizes is included, small power-law tem- 
perature dependences are predicted. Details such as the 
anisotropics (Ising-like or Heisenberg-like moments) as 
well as the mechanisms for dissipation are thought to de- 
termine the observability conditions for these effect o"'^^'^" . 
Fei_a;Coa;S2 appears to be a system where such effects 
may be observable over a wide temperature range as the 
Kondo temperatures we measure are small, of order 1 
K^^, while the Weiss temperatures are typically much 
larger, so that dissipation by way of the charge carriers 
may be relevant only at the lowest temperatures. 

It is no surprise that the disorder inherent to a doped 
semiconducting system can lead to the formation of an in- 
homogeneous magnetic state at low temperatures. What 
we find extraordinary are the unusual power-law temper- 
ature and magnetic field dependencies of physical prop- 



erties that we measure along with the very sharp, non- 
analytic, field dependence of the magnetic susceptibil- 
ity near H — 0. That a simple model of disordered 
magnets^^ which postulates the existence of rare ordered 
regions can describe the broad features of such a compli- 
cated system is both surprising and compelling. 



VII. ACKNOWLEDGMENTS 

We thank I. Vekhter and C. Capan for discus- 
sions. JFD, DPY, and JYC acknowledge support 
of the NSF under DMR084376, DMR0449022, and 
DMR0756281. We thank J. M. Honig for providing crys- 
tals of Vo.99Tio.01O3 for use as a manometer in our pres- 
sure cells. 



^ S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 

Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka- 

nova, & D. M. Treger, Science 294, 1488-1495 (2001). 
2 S. von Molnar & D. Read, Proc. IEEE 91, 715 (2003). 
^ H. Ohno, A. Shen, F. Matsukura, A. Oiwa, A. Endo, 

S. Katsumotos, & Y. lye, Appl. Phys. Lett. 69, 363-365 

(1996). 

^ H. Ohno, H. Munekata, T. Penney, S. von Molnar, & L.L. 

Chang, Phys. Rev. Lett. 68, 2664-2667 (1992). 
^ F. Matsukura, H. Ohno, A. Shen,& Y. Sugawara, Phys. 

Rev. B 57, R2037 (1998). 
® M. L. Reed, N. A. El-Masry, H. H. Stadelmaier, M. K. 

Ritums, M. J. Reed, C. A. Parker, J. C. Roberts, & S. M. 

Bedair, Appl. Phys. Lett. 79, 3473 (2001). 

R. C. Myers, B. L. Sheu, A. W. Jackson, A. C. Gossard, 

P. Schiffer, N. Samarth, & D. D. Awschalom, Phys. Rev. 

B 74, 155203 (2006). 
* G.Salis, Y. Kato, K. Ensslin, D. C. DriscoU, A. C. Gossard, 

& D. D. Awschalom, Nature 414, 6864 (2001). 
^ D. J. Priour & S. Das Sarma, Phys. Rev. Lett. 97, 127201 

(2006). 

^° V. M. Galitski,A. Kaminski, & S. Das Sarma, Phys. Rev. 

Lett. 92,177203 (2004). 
" T. C. Schuhhess, W. M. Temmerman, Z. Szotek, W. H. 

Butler, & G. M. Stocks, Nature Mat. 4, 838-844 (2005). 

S. Guo, D. P. Young, R. T. Macaluso, D. A. Browne, N. L. 

Henderson, J.Y. Chan, L.L. Henry, & J.F. DiTusa, Phys. 

Rev. Lett. 100, 017209 (2008). 
" H. S. Jarrett, W. H. Cloud, R.J. Bouchard, S. R. But- 
ler, C.G. Fredric, & J.L. Gillson, Phys. Rev. Lett. 21 617 

(1968). 

^'^ R. J. Bouchard, Mat. Res. Bull. 3, 563 (1968). 

^5 K. Ramesha, R. Seshadri, C. Ederer, T. He, & M. A. Sub- 

ramanian, Phys. Rev. B 70, 214409 (2004). 
" L. Wang, T. Y. Chen, C. L. Chien, J. G. Checkelsky, J. C. 

Eckert, E. D. Dahlberg, K. Umemoto, R. M. Wentzcovitch, 

& C. Leighton, Phys. Rev. B 73, 144402 (2006). 
" S. F. Cheng, G. T. Woods, K. Bussmann, I. I. Mazin, R. J. 

Soulen, E. E. Carpenter, B. N. Das, & P. Lubitz, J. Appl. 

Phys. 93, 6847 (2003); 

R. B. Griffiths, Phys. Rev. Lett. 23,17 (1969). 



E. Miranda & V. Dobrosavljevic, Rep. Prog. Phys. 68, 

2337-2408 (2005). 
^° T. Vojta, J. Phys. A:Math. Gen. 39, R143-R205 (2006). 

G. R. Stewart, Rev. Mod. Phys. 73,797-855 (2001) 

A. H. Castro Neto, G. Castilla, & B. A. Jones, Phys. Rev. 

Lett. 81,3531 (1998). 
2^ A. H. Castro Neto & B. A. Jones, Phys. Rev. B 62, 14975 

(2000). 

^'^ M. B. Salamon, P. Lin, & S. H. Chun, Phys. Rev. Lett. 
88, 197203 (2002). 

J. Deisenhofer, D. Braak, H.-A. Krug von Nidda, J. Hem- 
berger, , R.M. Eremina, V. A. Ivanshin, A. M. Balbashov, 
G. Jug, A. Loidl, T. Kimura, & Y. Tokura, Phys. Rev. 
Lett. 95, 257202 (2005). 

Y. Shimada, S. Miyasaka, R. Kumai, & Y. Tokura, Phys. 
Rev. B 73, 134424 (2006). 
^"^ J. Herrero-Albillos, L. M. Garcia, & F. Bartolome, J. Phys. 
Cond. Matt. 21, 216004 (2009). 

C. Ancona-Torres, D. M. Silevitch, G. Aeppli, & T. F. 
Rosenbaum, Phys. Rev. Lett. 101, 057201 (2008). 
S. Guo, D. P. Young, R. T. Macaluso, D. A. Browne, N. 
L. Henderson, J.Y. Chan, L.L. Henry, & J.F. DiTusa, sub- 
mitted to Phys. Rev. B (2009). 
^° L. Neel & R. Benoit, C. R. Hebd. Seances Acad. Sci. 237, 
444 (1953). 

^1 I. I. Mazin, Appl. Phys. Lett. 77, 3000 (2000); 

S. Guo, Ph.D thesis Mangetic, thermodynamic, and trans- 
port properties of the magnetic semiconductor Fei-xCoxS2 
and superconductor LaSb2. Louisiana State University 
(2006). 

K. Murata, H. Yoshino, H. O. Yadav, Y. Honda, & N. 
Shirakawa, Rev. Sci. lustrum. 68, 2490-2493 (1997). 
'^'^ K. Koyama, S. Hane, K. Kamishima, & T. Goto Rev. Sci. 
lustrum. 69, 3009-3014 (1998). 

K. Kamishima,M. Hagiwara, & H. Yoshida, Rev. Sci. lu- 
strum. 72, 1472-1476 (2001). 
'^^ S. A. Carter, T. F. Rosenbaum, M. Lu, H. M. Jaeger, P. 
Metcalf, J. M. Honig, & J. Spalek. Phys. Rev. B 49, 7898- 
7903 (1994). 

R.G. Goodrich, D. Hall, E. Palm, & T. Murphy, Cryogenics 
38, 221-225 (1998). 



22 



A. Arrott & J. E.Noakes, Phys. Rev. Lett. 19, 786 (1967). 
M.B. Maple, J. Wittig & K.S. Kim. Phys. Rev. Lett. 23, 
1375 (1969). 

J.F. DiTusa, K. Friemelt, E. Bucher, G. Aeppli & A.P. 

Ramirez, Phys. Rev. Lett. 78 2831 (1997). 

Y. Nishino, M. Kato, S. Asano, K. Soda, M. Hayasaki, & 

U. Mizutani. Phys. Rev. Lett. 79, 1909-1912 (1997). 

M. Lakner, H. v. Lohneysen, A. Langenfeld, & P. Wolfle, 

Phys. Rev. B 50, 17064 (1994). 

R. N. Bhatt and P. A. Lee, Phys. Rev. Lett. 48, 344 (1982). 
M. A. Paalanen, J. E. Graebner, R. N. Bhatt, & S. 
Sachdev, Phys. Rev. Lett. 61, 597 (1988). 
See e.g. S. Blundell, Magnetism in Condensed Matter Ox- 
ford University Press(Oxford) (2001). 

T. Westerkamp, M. Deppe, R. Kuchler, M. Brando, C. 
Geibel, P. Gegenwart, A.P. Pikul, & F. Steghch, Phys. 
Rev. Lett. 102 206404 (2009). 

See e.g. N. W. Ashcroft and N. D. Mermin Solid State 



Physics (1976) 

A. "j. Milhs, D. K. Morr, & J. Schmahan, Phys. Rev. B 66, 
174433 (2002). 

C. A. M. Mulder, A. J. van Duyneveldt, & J. A. Mydosh, 
Phys. Rev. B 23, 1384 (1981). 

K. Binderfe A. P. Young, Rev. Mod. Phys. 58, 801-975 

(1986). 

See e.g. Z. W. Ouyang, V. K. Pecharsky, K. A. Gschnei- 
dner, D. L. Schlagel, & T. A. Lograsso, Phys. Rev. B 74, 
094404 (2006). 

Wenhao Wu, B. EUman, T. F. Rosenbaum, G. Aeppli, & 

D. H. Reich, Phys. Rev. Lett. 67, 2076-2079 (1991). 

D. M. Silevitch, D. Bitko, J. Brooke, S. Ghosh, G. Aepph, 
& T. F. Rosenbaum, Nature (London) 448, 567 (2007). 
H. Hiraka & Y. Endoh, J. Phys. Soc. Jpn. 63, 4573-4582 
(1994). 



M (/Xb / Co) y (Q) 




o 



o 



o o 




o 



o 



8w". T, (K) 

h- ^ 

01 o 



I — ' • 

O 

o 

r-t- 



O 



o 
b 



o 
b 



> 
n 



o 




cn 



O 

O 
o 

'iX ^ 
cr cr cr cr cr cr S 

ffi p fD (I) (D 



1 



^ ^ ^ 



O 



o 
o 



o 
o 

f\3 



Curie Constant (K) 



o 



CJn 
o 
o 



cn 
O 
o 
o 



C/T (J/mole K') 



C (J/mole K) 




• )K X O O □ ♦ ; 

k> kvl k> s> k> • 

1^ k'N k'S 1^ 1^ . 

II II II II II II II 

o o o o o o o 

b o o b b b 

CX) CJi 4^ O O 



x<a» 



O 



cn 
O 







CO 

o 
o 



C/T (J/mole K') C/T (J/mole K') 



o 



o 



o 



o 



o o 



o 




o 
b 



• ° 


> 


+ 


♦ 


• X 


X 


X 


X 


• II 


II 


II 


II 


o 


o 


o 


o 



Ol 

Ol 



O O O O 
O O O 
fV) -s2 cn O 







■ 


II II 


X 


X 

II 


O O 


II 

o 


o 


b b 


b 


b 


CD CJn 






CJ1 





o 
b 



CD 

I — ^ 

>< ■ 

o 
o 

X 



o 



CO 



o 



o 
xb 



o 
b 

03 




2NJ/x (1/FU) 



S (xR ln(2)) 



o 



o 



o 




X 

o o 
b b 



♦ ► I 



O 



S (xR ln(2)) 




Wilson Ratio 



o o 



o 
o 



o 
o 

o 



► ♦ 

► ♦ 

► ♦ 

► ♦ 



♦ ■ 

► > 

► m 



o 

I — ' ' 

o 



CO 

X 



CD 
I — ' 

I 

X 

n 
o 

X 




CO 

CO 

o 



• 


■ 


♦ 


► 


X 

II 


X 

II 


II 


X 

II 


o 


o 


o 


o 


b 


b 


o 


b 






o 


o 






-n2 


cn 



o 



T, 



a 



o 



o 
b 



o 
b 



o 
b 



► ♦ o o ■ • 

Q P Q P P 
n 




M / (G^"0 



> A ^ o + 
CO -v2 ^ 02 
^ Oi cn 

□ • [> ► 

^ CD CO CD 



1^ I-^ 



M (G) 




□ B - 
II 

o °° 
cn ♦♦ 



135? -v2 en fo 

1^ ?N 



X A X V 

cn 4^ hji- 

i;s-< '-v2 cn Kd 
cn cn 



o 
b 



o 



o 

M 



o 



o o 



M (G) 



CO 



C/T 
(J/mole K') 



o 

O I 



O 

b 



O 

b 



a:) 



■ 


► 




X 

II 


X 


X 

II 


o 


o 


o 


o 


o 


b 




o 


o 






o 









C/T 
(J/mole K') 



o 
o 



o 
b 




I — ^ 
>< 

O 
O 

X 



X" (emu/mole) 



X' (emu/mole) 




O 

b 



o 



o 
b 



O 




o 



o 



1/x' 



o 



o 



o 



o 



X"{H)/x"(0) X'(H)/X'(0) 




